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Aims. We consider relativistic electrons in magnetospheric flows close to a supermassive black hole and study the mechanism of 

quasi-linear diffusion (QLD) to investigate the correlation between y-ray and radio emission in active galactic nuclei. 

Methods. Moving in the nonuniform magnetic field the particles experience a force that is responsible for the conservation of the 

adiabatic invariant. This force, together with the radiation reaction force, tends to decrease the pitch angles. Contrary to this, the QLD 

attempts to increase the pitch angle, and this maintains the synchrotron emission regime. To examine the balance between the QLD 

and the aforementioned dissipative factors we investigate the quasi-stationary state by applying the kinetic equation. 

Results. Considering the magnetospheric plasma close to the supermassive black hole, we examined the efficiency of the QLD for 

different parameters. By examining the cyclotron instability, we show that despite the efficient dissipative factors, the cyclotron modes 

excite transverse and longitudinal-transversal waves which leads to the QLD. We find that the QLD provides a connection of emission 

in the y-ray and radio domains. We show that under favourable conditions the radio emission from 22 MHz to 9 GHz on the mpc 

scale is associated with y-ray emission from 900 GeV down to 9 GeV on the same scale. 

Key words. Galaxies:active-lnstabilities- Magnetohydrodynamics (MHD)-Radiation mechanisms: non-thermal 



1. Introduction 

During the last decades the interest in very high energy (VHE) 
radiation of active galactic nuclei (AGNs) has considerably in- 
creased because of the detailed data from the new telescopes 
(MAGIC, HESS, EGRET, Fermi). Active galactic nuclei, and es- 
pecially blazars reveal a strong connection between high-energy 
y-ray and radio emission. Since the EGRET era, a direct con- 
nection between radio and y-ray flux density has been inten- 
sively discussed. In particular. Bloom (2008) performed a sta- 
tistical analysis of the broadband properties of EGRET blazars 
and showed a strong positive correlation between the total radio 
luminosity at 8.4 GHz and the y-ray luminosity in the EGRET 
band (Ly ~ LP/JJ). The comparison of the radio flux density at 
8.4 GHz and the y-ray (> 100 MeV) photon flux measured by 
the Fermi y-ray Space Telescope has been presented by Giroletti 
et al. (2010) and the Spearman's rank correlation coefficient was 
estimated to be 0.57. Kovalev et al. (2009) found a positive cor- 
relation between the pc-scale radio flux at 15 GHz and the y-ray 
photon flux in a sub-band of the Fermi/LAT detector (100 MeV- 
1 GeV), and a time separation of a few months between y-ray 
and radio flares. 

Most commonly the VHE emission of AGNs is 
explained in terms of the inverse Compton upscat- 
tering (Blandford et al. 1990) and curvature radiation 
(Gangadhara 1996, Thomas & Gangadhara 2005). The major 



reason for this is that if the magnetic field is strong, the syn- 
chrotron cooling timescale is very short and relativistic particles 
very soon transit to their ground Landau level, which prevents 
the subsequent emission process (IQsmanov & Machabeli 20lOl l. 



Recently, it has been shown that the quasi-linear diffu- 
sion (QLD) may create necessary conditions in the magne- 
tospheres of AGNs and pulsars for avoiding the synchrotron 
emission damping despite the strong magnetic field. According 
to the model of Kazbegi et al. (1992), the cyclotron insta- 
bility appears in the magnetospheres of pulsars. These unsta- 
ble modes lead to the feedback of the excited waves on par- 
ticles by means of the process of the diffusion. The pitch an- 
gles are no longer negligible and the emission process con- 
tinues (e.g. 'Machabeli & Usov 1979' 'Lominadze et al. 19791 
Malov & Machabeli 2001 Chkheidze & Machabeh 200 7). This 
model was applied to the Crab pulsar (Machab eli & 2009 I 
Machabeli & Osmanov 2010i) to explain the recent results from 
the MAGIC Cherenkov telescope between 2007 October and 
2008 February (lAIbert et al. 2002i 

In the context of works presented by Bloom (2008) and 
Giroletti et al. (2010) (where the observationally evident strong 
correlation of high-energy y-rays and low-frequency radiation is 
discussed), it is of fundamental importance that according to the 
QLD, two different radiation domains are excited. As was stud- 
ied by Osmanov & Machabeli (2010) and Osmanov (2010), the 
cyclotron resonance efficiently generates transverse momenta, 
which leads to non-zero pitch angles and the subsequent syn- 
chrotron radiation. This means that a connection between two 
radiation domains is naturally introduced into the quasi-linear 
diffusion process. In particular, examining the QLD in AGNs, 
Osmanov & Machabeli (2010) studied the generation of X-rays 
associated with the radio emission, and the work of Osmanov 
(2010) was related to high-energy emission associated with sub- 
millimeter/infrared radiation. 
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In this paper we study the production of high-energy radia- 
tion in the GeV-TeV interval in the magnetospheres of supermas- 
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sive black holes. The paper is organized as follows. In Section 
2 we present the model, in Sect. 3 we apply the mechanism to 
AGNs and in Sect. 4 we summarize our results. 



2. Theory 

In general, Lorentz factors of magnetospheric plasma particles 
m AGN^] lie in a broad interval ranging from ~ 1 up to ~ 10 
dOsmanov et al. 20 07 Rieg er & Aharonian 2008) , therefore, for 
simplicity one can consider the magnetosphere to be composed 
of the electron-positron plasma component with relatively low 
Lorentz factors, jp, and highly relativistic electrons, the so- 
called beam component electrons with the Lorentz factor ji, 

In this section we consider the QLD mechanism and derive 
an average value of the pitch angles driven by the diffusion. 
In general the process of energy dissipation that takes place in 
strong magnetic fields is dynamically controlled by two forces: 
one is the synchrotron radiative force (L andau & Lifshitz 1971| l 
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(1) 



where a - 2e^(y^/(3c^), is the pitch angle, llib = eB/mc 
is the cyclotron frequency, B is the magnetic induction, c is 
the speed of light, e and m are electron's charge and the rest 
mass, respectively; the other is the force G, responsible for 
the conservation of the adiabatic invariant, / = 3cp]_/2eB 
dLandau & Lifshitz 197 it 



mc 
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where p is the curvature radius of the magnetic field lines. In 
general, in the magnetospheres of AGNs two effects take place: 
the energy dissipation provoked by the forces, F, G, which lead 
to the decrease of the pitch angles, the QLD, which opposes the 
first effect. The dynamical process saturates when the effects of 
the above-mentioned forces are balanced by the diffusion, which 
leads to non-zero pitch angles and the subsequent radiation. 

As discussed by MachabeH & Usov (1979), Machabeli & 
Osmanov (2010), Osmanov & Machabeli (2010), for an activa- 
tion of the QLD, the cyclotron modes must be excited, which 
leads to the diffusion in terms of feedback of the cyclotron waves 
on particles. On the other hand, according to the work of Kazbegi 
et al. (1992), the anomalous Doppler effect can efficiently in- 
duce unstable cyclotron waves with the corresponding frequency 
dMachabeli & Usov 1979llOsmanov & Machabeli 20T0b 
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where 5 - cxPpliAuP'gY'p), ojp = ^]A■mpe^ Im is the plasma fre- 
quency and Up is the plasma density. Here we assume that energy 
in plasma is uniformly distributed, therefore, the plasma density 
can be approximated to be nbjhlyp, where ni, is the beam den- 
sity (Machabeli & 2009 ). Although the aforementioned expres- 
sion is a direct consequence of a resonance character of the cy- 
clotron instability, in exciting waves all resonance particles (with 
broad energy spectra) participate, therefore the range of spectral 
frequencies is wide and the waves are not characterized by a 



' The magnetosphere of AGNs with their typical lengthscales 10''*"''^ 
cm is a region around the supermassive black hole, where the magnetic 
field (lO'' G - close to the black hole and 10- 100 G - in the light cylinder 
zone) is dominant. 



monochromatic signature. In the framework of the present ap- 
proach, one has two simultaneously generated connected VHE 
y-rays and low-frequency cyclotron modes. This result seems to 
be promising in the context of recent investigations, because it 
shows the strong connection between the y-ray and radio fluxes 
dBloom 2 008, Gir oletti etal. 20101 1. 

A similar problem was studied by Osmanov & Machabeli 
(2010) who examined the physical regime \G±\ » \Fj_\ and 
|GJ <sc |FJ and studied the production of X-rays (0.13 keV-100 
keV) connected with the radio emission (0.04 MHz-35 MHz). A 
different regime, |G_l| <k \Fj_\ and IGJ «; |FJ, was examined by 
Osmanov (2010). This regime shows the possibility of produc- 
tion of MeV-GeV y-rays, which are strongly correlated with sub- 
millimeter/infrared radiation domains. In the present paper we 
examine the same regime, but for different parameters. Unlike 
the previous work, where the QLD was studied in the light cylin- 
der area where the magnetic field is of the order of ~ 300G, in 
the present paper the physics of the QLD is examined for a re- 
gion close to the supermassive black hole event horizon where 
B ~ IQ^G (see lThorne et al. 19881 ). Another difference is the fact 
that here we study the cyclotron waves generated by plasma par- 
ticles with Lorentz factors, jp ~ 2 - 4, whereas values of jp 
considered by Osmanov (2010) were higher by two orders of 
magnitude. As we will see, the aforementioned differences sig- 
nificantly change the results. By comparing the transverse and 
longitudinal components of the forces, one obtains 
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Evidently, the relations \G±\ <K \Fj_\ and IGJ <K IFJ are satisfied 
if i/r » 2.2 X 10"^ rad if we consider the physical parameters 
yt ~ 10\ B ~ 10^ G, Rg ~ p. 

By assuming a quasi-stationary scenario (d/dt = 0), the 
corresponding kinetic equation will be different from that of 
Machabeh & Osmanov (2010) 
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where / = /(i^, jb) is the distribution function of particles. 



(7) 



are the diffusion coefficients and \Ek\^ is the energy density per 
unit wavelength. 

For estimating \Ek^, we assume that half of the plasma en- 
ergy density, mc^nbybl2, converts to the energy density of the 
waves \Ek^k, then an expression of \Ek\^ writes as follows 
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Fig. 1. Behaviour of ejev with respect to ji,. The set of parame- 
ters is Mg = 1, B = 10"^ G, nt = 20 cm"^ and y,, e {2; 3; 4). 




Fig. 2. Behaviour of eTev versus v and jp. The set of parameters 
is Mg = 1 and B = 10"* G, rib = 20 cm^l 



For solving the aforementioned kinetic equation with respect to 
the pitch angles we separate the variables by expressing the dis- 
tribution function x{^)f{yb)- Then the equation governing the 
distribution by pitch angles is given by 



a^/ mcji, aip 



and the corresponding solution writes as 



where 



A = 



amcy^ 
AD,, ' 



(9) 
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(11) 



For simplicity we use the average value of the pitch angle, 
which after taking Eqs. (IIOII lb into account leads to 



J "CO 



Clearly, despite very efficient synchrotron losses the quasi- 
linear diffiision works against the dissipative forces and 
maintains non-zero pitch angles, which prevent the syn- 
chrotron mechanism from damping. Under these conditions, 
the relativistic particles with Lorentz factors, y, radiate in 
the synchrotron regime and emit photons with energies 



(Rybicki & Lightman 1979 1 



e,v ~ 1.2x lO'^Bj^sini//. 



(13) 



3. Results 



According to the observational evidence, AGNs are character- 
ized by emission in the VHE domain, which cannot be ex- 
plained in the framework of thermal processes. This can also 
be explained through the existence of ultra-relativistic parti- 
cles in the magnetosphere. In general, there are several accel- 
eration mechanisms that might provide very high Lorentz fac- 
tors. We show that processes such as Fermi-type acceleration 
jCatanese & Weekes 1999l l and acceleration owing the the black 
hole dynamo mechanism (ILevinson 20001 ) can provide Lorentz 



factors of the order of ~ 10^ close to the supermassive black 
holes, i.e. at r ~ Rg, where Rg = IGMbhIc^ ~ 3My x 10'"^ cm is 
the gravitational radius of the black hole, Mg = M/(10^Mo) and 
Mq is the solar mass. 

Osmanov & Machabeli (2010) have shown that the syn- 
chrotron cooling timescale for ultra-relativistic electrons in the 
magnetospheres of AGNs close to the black hole is of the order 
of ~ lO^^s-lO^^s. To estimate the efficiency of the energy losses 
it is reasonable to calculate the kinematic timescale fi,„ ~ Rglc. 
Considering the supermassive black hole with Mc, - 1, one can 
see that fA,„ ~ 10''s, which by many orders of magnitude ex- 
ceeds the dissipation timescale, therefore, the cooling process is 
extremely efficient and the electrons stop emitting very soon. 

Anomalous Doppler effect significantly changes the over- 
all radiative pattern of the system. Equation ^ clearly shows 
that the mentioned effect efficiently induces unstable cyclotron 
waves with the frequency (lOsmanov & Machabeh 2010] l 



(12) v.8.7xl0^(^) x(-) x(— ) x(-— )h.,(14) 



which in turn lead to the feedback of these modes to the particles 
and switch on the QLD, which generates the non-vanishing pitch 
angles and the subsequent synchrotron process. 

As a first example, we examine a supermassive black hole 
with M<) - 1 and the typical value of the magnetic induction for 
a nearby zone of the black hole B - lO'* G. Since we assume that 
energy in plasma is uniformly distributed, the value of plasma 
density ni,yi,/yp, for nb = 20 cm"^ and B - lO'* G is of the order 
of 10^ cm"^, which is compatible with observations of narrow 
and broad line regions and indicates that the plasma density in 
the environment of AGNs lies in the interval ~ (10^ - lO''*) cm"^ 
(IXu et al. 2007llArav et al. 2002l l. We consider the beam compo- 
nent's Lorentz factors to be of the order of 10** and the plasma 
component with the Lorentz factor, yp = 2. From Eq. (fT4l i one 
can show that the cyclotron frequency of the order of 87 MHz 
is excited, which in turn, by means of the QLD, leads to the av- 
erage value of the pitch angle ~ 0.2 rad (see Eqs. 01 11121 ) and 
the subsequent synchrotron emission in the VHE domain ~ 240 
GeV [seeEq. O]. 

In Fig. [1] we show the dependence of e^. on yb - The set of pa- 
rameters is Mc) = I, B ^ 10* G, rib = 20 cm"^ and yp e {2; 3; 4). 
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Fig. 3. Behaviour of exev versus v. The set of parameters is 
Mg = 1 and B = 10"^ G, rib = 20 cm-^ and y,, = {2; 3; 4). 
Evidentlythe curves are truncated, which is a direct result of 
the fact that the curves are one-dimensional sections of a two- 
dimensional surface. 




Fig. 4. Behaviour of erev versus y^. The set of parameters is 
Mg = 1 and B = lO'* G, y,, = 2 and tib = {0.2; 2; 20) cm^l 



Clearly, physically realistic parameters might guarantee the syn- 
chrotron emission from ~ 70 GeV (yt, = 10*^, yp = 4) to 0.9 
TeV iyh - 2 X 10^, yp = 2). These plots show a continuously 
increasing character of Cj.^^ with respect to yt,, which is a natural 
result of the fact that more energetic beam particles radiate pho- 
tons with higher energies. Unlike the beam Lorentz factors, the 
behaviour of e,,^^ with respect to yp is a continuously decreasing 
function. This can be seen from Eqs. (II 1112k ^ ~ ifDlI, which 
by combining with Eqs. ( 1718b reduces to ~ 1 /y^. Therefore, by 
considering higher values of jp, the resulting synchrotron emis- 
sion energy will be lower 

In Fig.|2]we show the behaviour of ejev with respect to v and 
yp. The set of parameters is Mg = 1 and B = lO'* G, = 20 
cm"^. The 2D surface shows that for lower values of the plasma 
Lorentz factors, the synchrotron emission energy is higher and 
the corresponding cyclotron frequency relatively is lower. Fig[3] 
shows the dependance of e^.^^ on the cyclotron radio frequency 
for three different values of plasma Lorentz factors yp = {2; 3; 4). 
Evidently, the VHE y-rays in the TeV domain are connected with 
the radio emission at ~ 87 MHz, and the higher frequency radio 
emission ~ 1.4 GHz is connected with the lower energy (~ 70 
GeV) y-ray synchrotron emission. 

One of the major magnetospheric parameters is the density 
of the beam component. In Fig.|4]we show the behaviour of e,,^^ 
versus n^. The set of parameters is Mg - 1 and B = 10** G, 
yp - 2 and - {0.2; 2; 20) cm"-*. Obviously, the higher the 
beam density, the higher the y-ray synchrotron radiation energy. 
As we already mentioned, the average value of the pitch angle 
is proportional to -^D^^, which by combining with D^j^ ~ 
[see Eqs. (17181 )1 confirms the results shown in the plots. For the 
considered interval of nt, the photon energy varies in the range, 
8.5 GeV-0.9 TeV. 

Fig Hpresents the dependance of e^^^ on v. From the plots we 
see that the TeV emission is connected to the radio band ~ 22 
MHz and the high-energy emission, 8.5 GeV is connected to 
~ 8.7 GHz radio domain. The present results are different from 
those of Osmanov & Machabeli (2010) and Osmanov (2010), 
where a different parametric space, has been considered. In par- 
ticular, Osmanov & Machabeli (2010) examined Mg - 1, y^ - 



2, yi, ~ 10^ ^, rth - 2000 cm ^ and Osmanov (2010) considered 
Mg = 1, yp e {200; 250; 300), yb ~ 10^ and «i = {5; 10; 15) 
cm"-'. 

According to the quasi-linear approach, the efficiency of the 
QLD depends on the efficiency of the cyclotron instability. For 
this purpose it is essential to consider the growth rate of the in- 
stability. Kazbegi et al. (1992) showed that the increment, F, of 
the modes is given by 



(ijj, 1 ui 

T^TT^ if 

ar/p 2 



(15) 



and 



2Lijyp 6 ■ 



1 ul 

if 

2 c2 
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where cDb = ^/47^^e^Jm is the plasma frequency of resonant 
(beam) electrons. For the parameters «/, = 20 cm"-', yp - 2 
(see Fig. [TJ, by assuming ~ c and p ~ Rg, one can show 
that M^./(2c^) <K 1, and the growth rate is defined by Eq. (fT5t . 
Therefore, the corresponding cyclotron timescale, fc,c = l/F, is 
of the order of 10 -' s, which is less by many orders of magnitude 
than the kinematic timescale, tki,, ~ Rg/c ~ 10^ s. The same rela- 
tion between the timescales is valid for the rest of the considered 
cases; therefore, the studied mechanism is extremely efficient. 

The interest in the QLD is twofold: on the one hand, 
the quasi-linear diffusion is an efficient mechanism that drives 
the synchrotron process despite very short synchrotron cooling 
timescales. On the other hand, as we have seen, the QLD pro- 
vides the possibility of generating connected VHE y-rays and 
radio emission on mpc scales. In particular. Bloom (2008) has 
found a high correlation of y-ray luminosity and radio luminos- 
ity. The work of Giroletti et. al (2010) is particularly interesting, 
because they show that radio emission at 8.4 GHz is connected to 
the VHE radiation (> 100 MeV). As we show in Fig.|5] the radio 
emission of the order of 8.4 GHz is associated with y-rays in the 
GeV band, therefore, the study of the QLD is very promising. 
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V (GHz) 

Fig. 5. Behaviour of ejev versus v. The set of parameters is Mg - 
1 and B - 10"* Q,jp - 2 and n\-, - {0.2; 2; 20) cm"-'. The curves 
are truncated for the same reason as in Fig. |3] 

4. Summary 

The main aspects of the present work can be summarized as fol- 
lows: 

1 . We studied the role of the quasi-linear diffusion with mag- 
netospheric plasma particles in AGNs. For this purpose we 
considered the kinetic equation that describes the balance be- 
tween the synchrotron reaction force and the diffusion. 

2. We examined our model for different values of physical pa- 
rameters. The QLD was studied versus the Lorentz factors 
of plasma and beam components respectively and versus the 
density of the beam component. We found that the higher 
the beam Lorentz factor, the higher the synchrotron emis- 
sion energy and moreover, that the higher the plasma Lorentz 
factor, the lower the radio frequency. We also showed that 
lower beam densities lead to lower energies in the y-ray syn- 
chrotron, and to higher energies in the radio domain. 

3. We also found that under favourable conditions the cyclotron 
instability generates radio emission from ~ 22 MHz to ~ 
9 GHz and the QLD provides synchrotron energies in the 
interval 900 down to 9 GeV. These results show the existence 
of connected radio and y-ray emissions on mpc scales. 



the acceleration timescale and the aforementioned instability 
timescale. We will investigate this problem in future studies. 
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The QLD is a mechanism able to explain a connection be- 
tween radio emission and VHE y-rays. On the other hand, a 
complete study requires an investigation of the spectral pattern 
provided by the mentioned mechanism. According to the stan- 
dard theory of the synchrotron emission fB ekefi & Barrett 19771 
[Ginzburg 1981 1, it is assumed that along the Une of sight the 
magnetic field is chaotic; correspondingly, the pitch angles vary 
in a broad interval (from to n) and the particles have uni- 
formly distributed pitch angles. Contrary to this scenario, in 
the framework of the present model the particles are nontriv- 
ially distributed (see Eq.fTOli. therefore, the spectral picture will 
differ from those of Bekefi & Barrett (1977), Ginzburg (1981). 
Another problem that has to be addressed is the variability of the 
corresponding VHE radiation. As the investigation shows, for 
the QLD to be efficient, two major requirements have to be ful- 
filled: the particles must be highly relativistic and the timescale 
of the cyclotron instability should be short enough. Therefore, 
the study of the variability implies the corresponding study of 



